Instability of strong regular reflection and 
counterexamples to the detachment criterion 

Volker Elling 

Abstract 

We consider a particular instance of reflection of shock waves in self- 
similar compressible flow. We prove that local self-similar regular reflec- 
tion (RR) cannot always be extended into a global flow. Therefore the 
detachment criterion is not universally correct. More precisely, consider 
the following angle condition: the tangent of the strong-type reflected 
shock meets the opposite wall at a sharp or right downstream side angle. 
In cases where the condition is violated and the weak-type reflected shock 
is transonic, we show that global RR does not exist. Combined with ear- 
lier work we have shown that none of the classical criteria for RR— >MR 
transition is universally correct. A new criterion is proposed. Moreover, 
we have shown that strong-type RR is unstable, in the sense that global 
RR cannot persist under perturbations to one side. This yields a definite 
answer to the weak-strong problem because earlier work shows stability of 
weak RR in the same sense. 

76H05; 76L05 



1 Introduction 

Consider compressible flow. In regular reflection (RR; see Figure [T|) an incident 
shock wave meets a solid wall in a reflection point and continues as a second, 
reflected shock. In many circumstances the effects of viscosity, heat conduction, 
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Figure 1: Left: local regular reflection (RR). Center: shock polar (t < t*, local 
RR possible). Right: r > t,, local RR impossible). 
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boundary roughness etc. are neglig so that inviscid models are appropriate; 
in this article we focus on compressible polytropic potential flow. Then shocks 
are sharp jumps satisfying Rankinc-Hugoniot conditions and the slip condition 
is used at walls: vi , V3 must be tangential. 

Consider a fixed constant state (velocity, density and sound speed) in the 

2- sector and vary the angle of the reflected shock. Each angle yields a different 

3- sector state. The curve of possible is called shock polar. The maximum 
angle between V2 and V3 is called critical angle. If it is larger than t (Figure 
[1] center), the angle between V2 and wall, then there are two possible reflected 
shocks satisfying the slip condition, called weak-type (W) and strong-typ^ (S). 

There is no local argument to rule out one type; the Rankine-Hugoniot 
and slip conditions allow both. At least in initial-value problems we expect 
uniqueness, in nature and in good mathematical models. For this we need to 
consider the global flow that contains the reflection, in particular domain shape 
and far- field /boundary conditions far from the reflection point. Of course there 
is an infinitco variety of such flows, but some observations and arguments apply 
to most if not all of them. 

If r > T*, on the other hand, then even locally RR is theoretically impossible 
because none of the reflected shock angles can make parallel to the wall 
(Figure[T]right). Around 1875, Ernst Mach [TU] discovered another pattern, now 
named Mach reflection (MR; see Figure IHl left), where incident and reflected 
shock meet off the wall in a triple point with a third shock, the Mach stem. 
For some parameters both RR and MR are possible. Starting with John von 
Neumann |llj . many researchers have tried to predict the precise parameters 
at which the RR^MR transition takes place (see [H [2] for a survey of this and 
other problems in shock reflection). 

There are three classical transition criteria. The von Neumann criterion does 
not apply in potential flo'wQ at all. The detachment criterion predicts global RR 
whenever a local RR exists. The sonic criterion, in contrast, predicts global RR 
if and only if there is a local RR with supersoni^ reflected shock. All three 
criteria are motivated by local considerations and well-defined for any global 
problem; of course the same criterion need not be correct for all global problems. 
However, we make a stronger observation: in a particular global problem, none 
of the classical criteria is correct, so that an entirely new criterion must be 
found. (The most promising candidates are modifications of the detachment 
criterion.) 

To define our problem wc add a second solid wall that meets the original wall 
right of the reflection point (Figure [2] left), enclosing a corner angle 180° — 9. To 
satisfy the slip condition in the constant-state 2-sector, the opposite wall has to 
move with horizontal speed w = w{9) so that V2 ■ n = w ■ n [n wall normal). 

^for experimental examples see 1141 p. 142f] 

^The names refer to their relative strength, but the absolute strength can be arbitrarily 
small or large. 

''in fact almost all flows with shocks include some form of shock reflection. 
*Even in Euler flow it applies only for sufficiently high Mach numbers, 
^which is necessarily weak-type 
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Figure 2: Left: trivial RR. Center: angle condition satisfied; may or may not ex- 
ist. Right: angle condition violated; no such flow can exist, even if the reflection- 
point tangent is allowed to be weak-type. 



There is exactly one 9 so that the new wall is perpendicular to the reflection- 
point tangent of the strong-type reflected shock (Figure [2] left). In this case, 
Vi ■ n = w ■ n as well, so the fluid in the 1-sector is also compatible with the 
wall. The result is what we call a trivial RR. 

However, for any other 6 the reflected shock would have to be curved (and 
border a non-constant region on its right), because its reflection point tangent 
does not form a right angle with the new wall. So there is a large variety of 
nontrivial cases; each has the same incident and reflected shock, but 6 and w 
vary. 

Alternatively, we may consider the coordinate system of an observer trav- 
elling in the wall-wall corner. He observes steady walls but moving shocks 
(Galilean invariance). Moreover, use reflection so that the new (opposite) and 
old (reflection) wall change places (Figure [3|) . 

Let a be the counterclockwise angle from incident shock to opposite wall in 
Figure [3] left. We have a family of problems, with parameter space consisting 
of triples {Mi, a, 9). At time t = the incident shock starts in the wall- wall 
corner (Figure [3] left). 

Such reflections occur in practice (Figure 3]). Experimentally a (nearly) 
straight vertical shock could be produced by breaking diaphragms or detonating 
small charges. This shock (Figure S] left) travels to the right through a tube, 
meeting the lower corner at some time. A first reflection occurs (Figure |4] 
left center). It is the classical case a = 0°, 6* < 90° which has been studied 
extensively [H [6l [HI [3]. The reflected shock travels up the wall, reaching a 
second corner at t = 0. In that instant, the local flow near the upper corner is 
the same as the initial data in Figure [3] left. 

[5] has already obtained global weak-type transonic RR for small perturba- 
tions of the trivial 9, in the following class: 

Definition 1. Consider self-similar potential flow (sec Section [2]). A transonic 
(or sonic) global RR (sec Figure [5] left) has a straight incident shock extending to 
infinity, meeting the reflected shock in a single reflection point on the reflection 
wall. The incident shock separates the 1- and 2-sector, two regions of constant 
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Figure 3: Left: initial data. Center: Weak-type trivial RR and possible per- 
turbations (dashed). Right: angle condition violated, neither weak-type nor 
strong- type global RR exist. 



Figure 4: Left: a = 90° incident approaching; left center: classical RR; right 
center: upper corner is locally like Figure [3l 



fluid state p, c, v. The reflected shock is C including the endpoints, separating 
the 2- from the 3-sector, meeting the opposite wall in a right angle. Flow in the 
interior of the 3-sector is elliptic (pseudo-Mach number L < 1, see (|10p ). with 
continuous fluid variables. 

The sonic criterion, in any reasonable precise formulation, predicts non- 
existence (and appearance of MR), so [Sj demonstrates that it cannot be uni- 
versally correct. The present paper considers the case of strong- type RR. 

Definition 2. We say Q satisfies the angle condition if the reflection-point 
tangent of the siroTi^-type reflected shock forms an angle < 90° (Figure [5] left 
and center) on its downstream side with the opposite wall. 

Theorem 1. Consider parameters AIi,a,6 so that local weak-type transonic 
RR exists, but the angle condition (Definition\^ is violated. Then global RR 
solutions of the kind in Definition]^ do not exist ( neither weak-type nor strong- 
type). 

Corollary 3. The following version of the detachment criterion is not univer- 
sally correct: 

^It is violated for a subset of the parameter space which is open and nonempty, hence 
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GenericalliP , when local RR exists, either weak- or strong-type can 
be extended into a global RR. 



Proof. It is sufficient to give a rigorous proof of existence of local RR satisfying 
the conditions of Theorem [TJ take x/t = to be the reflection point. Choose 
some supersonic 1-sector state. Then for sufficiently small r (Figure [T] left) we 
can find a weak- type incident shock and a 2-sector state with M2 > 1, as well as 
a strong-type reflected shock with V3 parallel to the wall. Choose an opposite 
wall whose extension to a line passes through the point x/t = V2 (so that the slip 
condition pT|) in the 2-scctor is satisfied) . If the angle between the two walls is 
chosen small enough, then the angle condition is violated. We can choose this 
local RR transonic as follows: by Proposition [5] applied to the incident shock 
polar, for sufficiently large r, i f , so r» J, for the corresponding reflected 
shock polar (Proposition [5]). For r « t*, M3 < f which necessarily happens as 
T grows. □ 

Corollary 4. Strong-type trivial RR is not always structurally stable. 

Proof. Given Mi, a, 6 for a strong- type trivial RR, perturb to a — d, 9 + 5 for 
some small 6 > 0. Then we arc in the situation of Theorem [1] where global RR 
cannot exist. □ 

In summary, we have obtained two separate results. First, the detachment 
criterion is not universally correct. Note however that wc have discussed only 
some cases with 9 > 90°. It would be interesting to find extensions to the 
classical case 9 < 90° = a; in that case, the detachment criterion is probably 
correct. A new RR— >MR transition criterion is proposed in Section [6l 

Second, while weak-type transonic trivial RR is structurally stable, .strong- 
type is not (Corollary |4]). This provides an important new answer to the weak- 
strong problem. Note that historically, dynamic stability, i.e. under perturba- 
tion of the initial data, has been considered. [7] observes numerically that both 
weak- and strong-type reflection are dynamically stable, so any mathematical 
result to the contrary appears to use an overly restrictive definition of stability. 



2 Potential flow 



Self-similar potential flow is the second-order quasilinear PDE 



V • (pVx) + 2p = 0. 



(1) 



Here 



(2) 



"generic" by any reasonable definition. 
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^ = 77) = x/t are similarity coordinates, x is called pseudo-potential. ip is the 
velocity potential: physical velocity is 

V = VtA- 

Moreover, density is 

P = ^-'(-X-^|Vxn. (3) 



TT satisfies 



where 



|=,-.|..-V^ c- (4) 



is the equation of state {po, Cq free parameters). The ratio of heat^ is restricted 
to 7 £ (1, 00). Differentiation of ^ yields the non-divergence form 

(c2/ - - C?) : = 0. (5) 

Here A : i? is the Frobenius product tr(A-'"i?), w(E)w ^ ww'^ (as opposed 

to \w\'^ ^ w ■ w) and is accordingly the Hessian. In coordinates: 

(c2 _ (^^ _ ^)2)^^^ _ 2(^^ _ ^)(^^ - r^)^^^ + (c2 _ (^,^ - r;)2)V',, = 0. (6) 
c is the sound speed, defined by 

c2 = c^ + (l-7)(x + ^|VxP). (7) 
It is sometimes more convenient to use the form 

{c^I - Vx') : V\ + 2c' - |Vx|' = 0. (8) 

This form is manifestly translation-invariant. Translation is nontrivial: in 
{t, X, y) coordinates it corresponds to a change of inertial frame 

v^v — w, ^ = x/t^^ — w, (9) 

where w is the velocity of the new frame relative to the old one. Obviously the 
pseudo-velocity 

z ;= Vx = V?A - ^ 

docs not change. 



^also: isentropic coefficient 
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Self-similar potential flow is mixed-type; the local type is determined by the 
coefficient matrix c^/ — Vx^ which is positive definite if and only if L < 1, where 

L := M ^ (10) 

c c 

is called pseudo-Mach number. For L > 1 the equation is hyperbolic; parabolic 
is L = 1 . L and z are the Mach number and velocity perceived by an observer 
traveling on the ray x = tS^. 

On a solid wall the slip condition 

Vx-n = (11) 

holds; for an observer traveling on the wall it corresponds to the usual 

■y. n = V-0 ■ n = (12) 

3 Shock conditions 

The weak solutions of potential flow are defined by ([T]). The corresponding 
Rankine-Hugoniot condition is 

PuZu = Pdz2 (13) 

where u,d indicate the limits on the upstream and downstream side and z", z* 
are the normal and tangential component of z. As the equation is second-order, 
we must additionally require continuity of the potential: 

= i^d- (14) 

By taking a tangential derivative, we obtain 

4 = 4 =: (15) 

Observing that cr = ^ • n is the shock speed, we obtain the more familiar 
form 

PuVl - PdVd = o{Pu ~ Pd), (16) 

(17) 

Fix the unit shock normal n so that z^' > which implies > as well. 
To avoid expansion shocks we must require the admissibility condition z^ > z^, 
which is equivalent to 

vl > vl- (18) 

We choose the unit tangent t to be 90° counterclockwise from n. 

By (fT7|) the tangential components of the velocity are continuous across 
the shock, so the velocity jump is normal. Assuming > (positive shock 
strength), we can express the shock normal as 

(19) 

\Vu - Vd\ 
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Figure 5: Left: transonic global RR, angle condition satisfied. Self-similar po- 
tential flow is elliptic in 17, hyperbolic elsewhere. Right: Shock polar argument. 

4 Nonexistence of some global RR 

We start with some facts about the shock polar. 

Proposition 5. Consider arbitrary Cu, Pu > and Mu G (l,c>o) and set Vu = 
(M„c„,0). For each (3 S (—90°, 90°) there is a steady shock with downstream 
unit normal fi = (cos /3, sin Its downstream data depends smoothly on /3. Let 
T he counterclockwise angle from to Vd- We restrict \(3\ < arccos -p- so that 
the shock is admissible. 

Then the shock polar f3 ^ Vd is smooth and strictly convex, with dpVd 
nowhere zero. 

There is an angle r* £ (0°,90°) so that each r G (— t*,Th.) is attained for 
two different (3. The one with smaller \vd\ yields a strong-type shock, the other 
one weak-type. For \t\ = t* they are identical and critical-type. 

There is a Ts & (0,t,) so that the weak-type shocks are supersonic for \t\ > 
Ts, transonic for \t\ < Tg. The other types are always transonic. 

If Mu I 1 with Pu, Cu fixed, then | 0. 

Proof. Most has been shown in |51 Theorem 1] and Proposition 2.10]; we 
only need to prove the last statement. Admissible shocks are those for \(3\ < 
arccos-j|-. As M„ J. 1, this range shrinks to {0}. By continuity, all points on 
the shock polar approach In particular t» j 0. □ 

Let f2 be the 3-sector excluding boundary (Figure [5] left). A opposite wall, 
B reflection wall. S reflected shock, each not containing its endpoints. 

Proposition 6. Consider the setting of Theorem [II with (0, 0) the wall-wall 
corner (see Figure\^. The vertical straight shock with upstream data V2,P2,C2 
through the reflection point has a downstream velocity wo = (vq,0) with 
iiq > 0. The same holds for all vertical shocks to the right of it. 
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Proof. First change to a coordinate frame with origin in the reflection point. 
For this observer the 2-sector velocity is V2 — £,r, which points into the reflection 
wall (FigurcOright). Consider shocks with that upstream velocity and upstream 
density p2 and sound speed C2; let Vd ~ be the downstream velocity. Let /3 
be the counterclockwise angle from V2 — S,r to shock downstream normal n; by 
(|19p fi is a positive multiple of V2 — Vd- 

By design (slip condition), the velocities — and Vg — for weak- type 
and strong-type reflected shock arc on the extension of the reflection wall into 
a line (sec Figure O right). By assumption of Theorem [1] the angle condition 
is violated, so the strong-type and therefore the weak-type reflected shock tan- 
gent in the reflection point arc down and strictly right. Thus the vertical shock 
through the reflection point has smaller \f3\ than either type, so by strict con- 
vexity of the shock polar (Proposition [5]) vq — ^r points into the reflection wall 
(see Figure [5]). Therefore Vq > 0, since V2 and thus vq arc horizontal. 

By [SI Proposition 2.9], vertical shocks more to the right have — ^r > 
Vq — £,R (because they are weaker, so Vd is closer to W2). Hence > Vq > as 
weU. □ 

Proof of Theorem [H Consider the same coordinates as in the statement of Propo- 
sition [51 Restrict ip to fl (taking its fi-side limits on dU). Let ipo be the value 
of ■0 in the reflection point ^n. Let Sq be the straight vertical shock through 
the reflection point; let ctq be its ^ coordinate. 
Consider a transonic global RR. 

Again by assumption the angle condition does not hold, so the reflection 
point shock tangent points down and strictly right (as in Figure [3| right, as 
opposed to Figure [5] left) . The upstream velocity Vip = V2 has ipx = V2 > 0, 
so necessarily -0 > -i/^o at the shock near the reflection point S^r. Therefore, the 
global maximum of ip over 57 (which must be attained since $7 is compact and 
ip continuous) is greater than ipQ and not attained in £^r. 

Consider a maximum > f/^o in a point ^ G S" — {^_r}. The shock tangent 
is vertical in ^ (by -04 = for a maximum at 5*; by the slip condition (|12p at 
A for a maximum in the point where S meets A). Moreover, tp > ipo implies 
the shock is right of the vertical reflection point shock because 1/' is continuous 
across the shock and ipx = V2 > on the upstream side. Hence by Proposition 
^ tp^ = > Vq > 0. This is incompatible with a local maximum. 

Hence ip does not attain its fi-maximum anywhere on S. Then the same is 
true for 

4>:=i' + 6^ 

if we choose S > sufficiently small. By linearity 

(/ - c-^'^x^) ■■ v^^/' = (/ - cr^Vx^) ■■ = 0, 

so by the strong maximum principle the maximum is not attained in Q either 
(if we choose J > so small that tp, like tp, cannot be constant). Moreover 

\/ip-n = \7ip-n + Srf =5rf >Q 
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Figure 6: [Reflection and opposite wall are exchanged in these diagrams.] Left: 
MR for 9 = 147.9°. Instead of meeting on the bottom wall, incident (left) and 
reflected (right) shock meet in a triple point with a vertical shock (Mach stem). 
Right: RR for 9 =■ 137.9°. Near the reflection point the flow is hyperbolic; the 
transition to elliptic is discontinuous (MR). 

on A and B, so the Hopf lemma rules out local maxima there. 

Finally, the boundary conditions on A and combined with continuity 
in (Definition [1]), imply VV'(O) = 0, so ?/'^(0) = V^(0) + 5 > 0, thus a local 
maximum in is impossible. 

We have ruled out every possible global maximum point in il. The contra- 
diction demonstrates that no with the desired properties exists. □ 

5 Numerical comparison 

Theorem [T] concerns the range of parameters with transonic weak-type RR, 
which is so narrow (see Figure [7] right) that numerics and experiments have not 
been able to settle questions for these flows. However, the range with super- 
sonic weak-type RR violating the angle condition is much larger and certainly 
interesting by itself. 

For 7 = 7/5, Mi « 3 and a = 0°, = 142.9° corresponds exactly to a 
strong- type trivial RR (i.e. strong- type shock perpendicular to opposite wall). 
We change by 5° to 147.9° without changing -\- a ot Mj. This way the 
opposite wall angle changes, but not the local RR parameters. The numerical 
results in Figure [6] left show an MR. 

We also study the opposite perturbation, to 6' = 137.9° (see Figure [6] right) . 
As expected there is still a local RR. The shock is essentially the strong-type 
shock, except in a small neighbourhood of the reflection point where it is weak- 
type and slightly hyperbolic. As t 142.9°, this neighbourhood shrinks to 
zero; it appears that the pattern converges to the trivial strong-type RR in 
this manner. This is why strong- type reflections are observed at a large scale 
sometimes. Note that there is a MR as well: at the transition from hyperbolic 
to elliptic. 

The calculations were made with a second-order scheme on an unstructured 
grid; other choices have no influence on the qualitative structure (RR vs. MR). 
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Figure 7: Proposed RR^MR transition, for 7 = 7/5 potential flow with a = 
0°. Right: detail. Weak-type reflection is supersonic above the "sonic" curve, 
transonic below; neither type exists below the "detach" curve. Theorem [T] rules 
out global RR below the solid curve. 



In principle Definition [T] and Theorem [T] could be extended to the supersonic 
cases. But while in the transonic case the flow is simple and predictable, at 
least for small perturbations from trivial RR, the supersonic cases can have 
several different qualitative structures. [U [6j construct self-similar RR with a 
continuous transition from hyperbolic to elliptic in the 3-sector, but Figure [6] 
right shows a MR, i.e. a discontinuous transition; double Mach reflection and 
other more complicated flows are possible too. Proving nonexistence in function 
classes large enough to accomodate all these structures is far beyond present-day 
techniques. 

6 Interpretation 

Despite the theorem and numerical examples, it is likely that the detachment 
criterion is still valid over a large part of the parameter space. In particular, 
the author believes that it is correct in the classical case a = 90°, 9 < 90°. We 
propose the following new criterion: 

The global flow is RR if and only if local RR exists and angle 

condition is satisflcd. 
Strong-type RR would appear only in the trivial right-angle borderline case 

separating global RR and global MR. 

Figure[7]shows the regions predicted by this criterion for 7 = 7/5 and a = 0°. 

Cases that have already been treated by construction of an exact solution or 

another rigorous method: 

1. Nonexistence of global RR below the detachment criterion is trivial. 

2. Nonexistence of global RR below the solid curve is done in this article for 
transonic weak- type RR by Theorem [T] 
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3. Existence of global transonic RR is done by [5] for some neighbourhood 
of each point on the transonic part of the "weak trivial" curve in Figure 
[7] (excluding endpoints) . 

4. [U [S] construct global supersonic RR for some of the supersonic param- 
eters, in particular some neighbourhood of each point of the supersonic 
part of the "weak trivial" curve in Figure [7] (excluding endpoints). 

In principle, [H [51 [HIS] go a long way towards constructing global RR in all cases 
not covered so far. In comparison, global MR is very difficult: the triple point, 
well-known to be theoretically impossible (von Neumann paradox), has a very 
complicated detail structure, according to numerical results of Hunter/Tcsdall 
(BM), see also [Hill]), 
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